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Abstract 

We propose a multidomain spectral collocation scheme for the approximation of 
the two-dimensional Stokes problem. We show that the discrete velocity vector field 
is exactly divergence-free and we prove error estimates both for the velocity and the 
pressure. 
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in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 


l 





Introduction. 

Domain decomposition methods for the approximation of partial differential 
equations are based on the partition of the physical domain into subdomains. These 
techniques allow one to decompose the problem into a collection of problems of smaller 
size, defined on the subdomains. By using proper interface conditions the overall ef- 
ficiency of the algorithm 'is then improved. 

If spectral schemes are considered domain decomposition methods may also be 
used to handle problems defined on plurirectangular domains. Actually, spectral 
methods are naturally defined only on rectangular domains (see e.g. [13]). If a 
problem on a plurirectangle is considered then it is possible to decompose such a 
domain in rectangles and employ spectral schemes on these subdomains. For the 
analysis of spectral multidomain methods for the Helmholtz equation we refer to 
[7], [9], [ll], [12], [16], [18], [19]. General considerations on domain decomposition in the 
framework of spectral methods can be found in [8, Ch.13]. 

In this paper we consider the Stokes problem in a plurirectangular domain with 
homogeneous Dirichlet boundary conditions. We analyze a scheme obtained from a 
generalization of the spectral collocation method introduced in [6] and [15] for the 
case of a monodomain approximation. An interesting feature of the numerical scheme 
is that, for any value of the discretization parameter, the approximated velocity field 
is exactly divergence-free. This is an important property in itself because it ensures 
that, even if a small number of degrees of freedom is used, the continuity equation 
is solved exactly. Furthmore, if the Stokes solver is used at each step of a time 
discretization of the full Navier-Stokes equation, then this property makes the time- 
algorithm more stable. 

Our main interest lies in the analysis of the discrete problem and in the proof of 
error estimates both for the velocity and the pressure. The number of subdomains is 
considered fixed and convergence is achieved when the degree of the solution (which 
is a polynomial in each subdomain) tends to infinity. Iterative methods for solving 
the discrete problem are under investigation and some related numerical experiments 
are being conducted. 

A multidomain scheme for the approximation of the Stokes problem can be found 
in [16] where the "Spectral Element” method is presented. The main difference be- 
tween that method and the one we analyze here lies in the choice of the discrete 
pressure space and in the quadrature formulae which axe employed. Domain de- 
composition for solving the Stokes problem in the framework of the Finite Element 
method is treated e.g. in [17] and [14]. 

The paper is organized as follows. In section one we introduce the continuous 
Stokes problem. A first version of the discrete problem is then formulated in a 
variational form and interpreted in a collocation one. Throughout the section we 
assume, for the sake of simplicity, that the physical domain is a rectangle partitioned 
in aligned subrectangles. We show that the discrete velocity is exactly divergence-free 
and we state an inf-sup condition which garantees compatibility between the velocity 
space and the pressure space (see [2]). In section two we give a new formulation of 
the discrete problem, which is more appropiate for the computation of the pressure. 
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Using results from the previous section, we prove convergence results both for the 
velocity and the pressure. In section three we consider the case of a plurirectangular 
domain and we extend the results of the previous sections. Finally an appendix is 
devoted to the proof of the inf-sup condition stated in section one. 


Notation. 

In the following M is a fixed positive integer and Cl is a plurirectangle partitioned 
in rectangles flj, i = 1, ..., M. Any quantity defined in fl» is identified by a superscript 
index i. For instance if v is a function defined on Cl then the restriction of v to 
Cli is denoted by u*. We shall use the classical Sobolev spaces H* (fl), H*(Cli), 
i = 1 s € R (see e.g. [l]). The norms of these spaces will be denoted by 

|| • ||, and || • ||,, i respectively. 

For any positive integer N we denote by Pm (fl,) the space of the restrictions to 
Cli of the polynomials of degree less or equal to N with respect to each variable. For 
any integer t, 1 < i < M, let 5^ be the set of the ( N + l) 2 nodes related to the 
Legendre Gauss Lobatto quadrature formula in Cli. The points of E X N are denoted by 
(x£,yj), 1 < kjl < N and the related quadrature formula is 


N N 


» 1 * IV 

/ v{x,y)dxdy a EE u (4>y/)<4“/> 

Jcii k=0 1=0 


( 0 . 1 ) 


where and uj] axe positive weights (see e.g. [10]). We recall that if v 6 P 2 Ar-i(n») 
then the quadrature formula (0.1) is exact. For any v,w € C°(n^) we set 

N N 

(v,w)N,i = 52 S w (*fc»»f) w fc w i» * = !»•••» Af- (0.2) 

k=o 1=0 

This bilinear form is an inner product on Pm (Cli), moreover we have (see [8]): 

Gi||v||o,» < ||v|jjv,< < C 2 ||t;||o,i Vv € P^(Ut), * = 1,..., M, (0.3) 

where C i and C 2 are positive comstants independent of N and || • ||jv,; is the norm 
associated to (0.2). For any v,w E C°(Cl) we define 

M 

({v,w)) N = ^(uSu/^j v,i. (0.4) 

*=1 

This bilinear form is an inner product on -Ptf(Di). Denoting by || • [|jv the 

associated norm, we have 

M 

Ci\ |v||o < \\v\\ N < C 2 \\v\\o Vv € (0.5) 

»=i 
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where Ci and C 2 are positive constants independent of N, M and of the measure of 
each subdomain. If v and w are vector valued functions we keep the same notation 
for the inner product and the norm induced by (0.2) and (0.4). Note that the form 
((•,•))// is an approximation of the L 2 -inner product denoted by (•,•)• 

Given an integer »\ 1 < i < M, and a point P_= (4,y?) € the characteristic 
polynomial associated to P is the element of Pjf (fl») such that 


Xp = 


in P 

in the other points of E^. 


( 0 . 6 ) 


Finally if is a rectangle given by x [y<— i,y*] we define, for any n e N: 

(0.7) 


^(y) ^ 


where L n {£) is the Legendre polynomial of degree n defined in [-1,1], (see e.g. [10]). 
We also let C denote a positive constant. The value of C shall always be independent 
of the discretization parameter N . 
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The continuous and the discrete problems. 

In this section we consider the case of a rectangle fl partitioned in aligned rect- 
angle fl*, i = 1 , M, as in the figure 1.1. The interfaces between subdomains are 
denoted by I 1 *, * = l,...,Af — 1. 



figure. 1.1 Decomposition with aligned subdomains. 


Let us consider the Stokes equations: 


— Au + Vp = f in fl 


V • u = 0 in fl (l.l) 

u = 0 on dfl. 

As usual u and p denote the velocity and the pressure of an incompressible fluid, 
the kinematic viscosity has been normalized to one and f represents a source term. 
Supposing f € (Ir 2 (n )) 2 it is well known that problem (l.l) has a unique solution in 
(^(n)) 2 x and that its variational formulation is the following (see e.g. 

( 21 |)- 

Find u € (j?o(fl)) 2 and p 6 L 2 (fl)/R such that 


( Vu, Vv) - (p, V • v) = (f , v) V v € (ff£ (fl)) 2 

(V -u,q) = 0 V 9 <E L 2 (fl)/R.- 


( 1 . 2 ) 


We now define a multidomain spectral approximation of this problem. The ve- 
locity space (ifo(fl)) 2 is approximated by its finite dimensional subspace V N defined 
as follows: 

V N = {v 6 (C°(fj)) 2 : v* e{P N (Cli)f, v* =0 on dn ndtli, i = 1 , ...,Af|. (1.2) 

Denoting by Mn the discrete pressure space (which is defined below) and recalling 
the definition (0.4), we consider the following approximation of problem (1-2). 
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(1.4) 


Find un € Vs and pn € Mn such that 

((Vuat, Vvjv))^ - {{pn, V • vn))n = ((f,vjv))Ar V w N € V N 


((V • u^r,gjv))jv — 0 

Here we have supposed that f € (C°(n)) 2 . We require that Mn is a subspace of 
n£iftr(fc ) such that problem (1.4) is well posed. Clearly Mn must not contain 
any “spurious mode” of the pressure, i.e. any nonvanishing element of the space 


M 

Zn,M = U e Jl P N{^i) : ((?>V' v ))w=° Vv<=V W }. 

»=1 

In the following ?pr position, which is proved in the appendix, we precise which 
are the elements of Zn,m- Let us first remark that, due to the geometry of Q, the 
polynomials l4(y), defined in (0.6), do not depend on t; hence the index * will be 
omitted. 

Proposition 1.1. The space Zn,m is a linear subspace of Il£li^(^») w ^ ose 
dimension is 4M + 4. A basis of Zn,m is given by the following elements. 


(xio.O 0), (xi,,o,...,o) 


(0 0,x2L), (°»-»o, X&) 


( 1 . 6 ) 

(1.7) 

(0, ..., 0, Xa ; > — XaT l > •••» ®)» (0, ..., 0, Xbi > — Xbt > 0, 0) for * = 1, ..., AT — 1 (1*8) 

(1,1,..., 1), ( L N (y),L N {y),...,L N (y )) (1.9) 

(0, ..., 0 , L* n (x) , 0 , ..., 0) , (0,...,0,L^(x)L w (y),0,...,0) for * = 1, ...,M. (1.10) 


We shall define the space Mn so that Mn © Znm = IIt=x Pn{IU)* a possible 
choice is 

M 

Mn = S JJ Pn{^x) : ((9.»‘))w ==:0 Vr € Zn,mJ- (l* 11 ) 

i=i 

Using (0.2) (0.4) and (l.6)-(1.10) it is readly seen that Mn can be written as follows: 


M 

M n = {? € II M&) : q 1 {ao) = Q l {*> o) = = 7 M (M = 0, (1.12) 


i=i 


- 1 , 

M M 

5^(9*. i )n,* = y = o, 

1 i=i 

(7*> L*n ( a: ))r/,i = 0, (7*, &n {x)Ln { y))N,% = * = 
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Theorem 1.1. Problem (1.4) has a unique solution. 

Proof. If f = 0 then choosing Vat = Uj\r and = p n in (1.4) we obtain 
((Vutf, Vujy))jv = 0. From (1.3) and the Poincar4 inequality this gives U// = 0. 
From the proposition 1.1 we also deduce pat = 0 hence the theorem is proved. I 

For the couple of spaces (Vat, Mat) the following “inf-sup” condition holds (see 
the appendix for the proof). 

Theorem 1.2. There exists a positive constant C independent of N such that 


• f ((gAT, V • Vat))aT > C_ 

q»eM N Vw€ v w IIgaHIat l|^ v Ar||iV — N 


(1.13) 


I 

Remark 1.2. Note that the right hand side of (1.13) depends on the discretization 
parameter N. This estimate is optim al, actually an upper bound for the left hand 
side of (1.13) is Cl/N , where C7 is a positive constant independent of N. For a 
detailed discussion on this subject we refer to [20]. | 

Proposition 1.2. If v^y € then V • vat G Mat- 


Proof. Since vat = 0 on dU we have V • v = 0 at oo, bo, om , and b^. Moreover 
v* = v* +1 on Ti for * = 1,...,M — 1, hence V-v*^) = V-v* +1 (a»+i) and V-v*(6i) = 
V • v* +1 (6j), for i = 1, ..., Jlf — 1. Recalling (0.2) and the fact that the Legendre 
Gauss Lobatto formula is exact for polynomials of degree 2 N — 1, we have, from the 
divergence theorem 

M 

E(v-v‘iw=°, 

ial 


moreover 


M M 

• v\L N {y)) N<i = ^K iie ,LAr(y))Ar,i - 


1=1 


1=1 


M N 


™ fa* 

= L n{.Vi) / v[ tX {x, yi)dxu\ = 0. 

»=1 1=0 


Using the same arguments we also deduce 


* rbi 

* V j Lff (x) ) AT,t' = ( u 2 ,y i {x)) N ,i = / „ I'M (^fc ) / u 2,y (®fc» y)dywj. =0 

k=0 


Finally, since v‘ e (P^iCli)) 2 we have (V • , L i N (x)L N (y)) i =0. It is readly seen 

that this implies (V • (y ))//■,» = 0, for * = 1,...,M. 

I 

The following is a direct consequence of (1.4) and Proposition 1.2. 


6 



Theorem l.S. If vln is the solution of problem (1.4) then V • Ujy = 0 in ft* for 

* = I 

Remark l.S. Using the technique introduced in [6, Corollary 5.1], for the mon- 
odomain case, one can prove that Ms is exactly the image of V N through the diver- 
gence operator. I 

We now give a collocation form of problem (1.4). Let us recall that B.* N n flj, 
for 1 < * < M is the set of the collocation points internal to fl», and E l N n I\, 1 < 
t < M — 1 is the set of the points which lies on the interfaces (not including the 
endpoints a* and 6»). Let now ££, k = 1,...,4, be four points of n Hi and let 
k = 1,2 be two points of 5^ n fy, for i = 2 We suppose that 

detjv?, 1 ^ )} ^ 0 (1.14) 

where <pj, l = 1, ..., 4 are the four polynomials 1, L l N (x) , L N (y ) , L l N {x)L^{y). More- 
over denoting by tf)j, l = 1,2, i = 2 ,...,M the polynomials L' N [x) and L x N (x)L N (y), 
we suppose that 

det{#’(f fc )} # 0 for i = 2 (1.15) 

The hypothesis (1.14) and (1.15) are for instance verified if for any i = 1 the 

points Cf, l = 1,2 (or / = 1,...,4 if i = l), axe not aligned. 

We can now state the following. 

Proposition l.S. Problem (1.4) is equivalent to find u n £ V) v and pu E Mjj 
such that 

-Au^(0 + Vp^(0 = f{C) ve € S/y n Cli for i = 1, ...,M, (1.16), 

V • u]v(£) = 0 V£ € 2 l N exept C € {a 0 ,&o,«iA*£fc> k = l,-> 4 }, (1.17) 

V-u}y(£)=0 E Sjy exept C € {<*»»&»> & = 1> 2 } for * = 2,..., M, (1.18) 

= ( 1 - 19 ) 

= - { (- Auj, + Vp*, - f>J, } U) Vf € Hj, n I\ < = 1 M - 1. 


■ 

We omit the proof of this proposition because it is a direct generalization of the 
proof given in [6, Prop. 5.1] for the monodomain case. Let us however precise some 
remarks about (1.16)- (1.19). In (1-16) we impose the momentum equation at the 
internal collocation nodes of each subdomain. These equations are obtained by (1.4) 
by choosing v^r as the characteristic function of the internal nodes of 0* , * = 1, ..., M 
and by using (0.2). In (1.17) and (1.18) the continuity equation is imposed in all the 
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collocation points exept a set of 4 M + 4 points. In order to have the same number 
of unknowns and equations it is necessary to eliminate such points (we recall that 
M*r is a subspace of codimension 4 M + 4 of flili. *V(n» jj- We enphasize that the 
solution of (1.16)-(1.19) do not depend on the choice of the points ££, k = 1,...,4 
and (? k , k = 1,2, * = 2 ,...,M, provided the hypothesis (1.14) and (1.15) hold. 
Concerning (1.19), note that on the right hand side we have the residual of the 
momentum equation multiplied by the weight w J N (which behaves like N~ 2 ). Hence 
(1.19) asymptotically ensures that the stress tensor is continuous at the interfaces 
T if i = l, ..., M — 1. This property is clearly verified by the solution of the continuous 
problem (1.1). 

2. Error estimates and a new pressure space. 

In this section we first prove an error estimate for the velocity field. 

Let w € (.Ho(n)) 2 , (where fl is given in figure 1.1), be such that V • w = 0 in 
H. We denote by n^ d * w w the orthogonal projection of w, for the inner product of 
(i?o(n)) 2 , on the space of the divergence-free functions which belong to V^. The 
following approximation result holds. 

Theorem 2.1. Let w € (Hq (n)) 2 be a divergence-free function and suppose that 
w* € (#'• (Hi)) 2 , with Oi > 1, for » = 1 There exists a positive constant C 

independent of both w and N such that 

l|w - nj»'*'w || 1 < cf; N l ~ ai ||w‘ !!„,(. (2.1) 

t=i 

Proof. Since w € {Hq (D)) 2 and V • w = 0 in fl there exists 0 e -ff^(fl) such 
that w = rot 0. In [5, thm. 3.6] it is proven that there exists 0// € 3$(Q) such that 

0V € 1 on Ti and on for * = 1,...,M - 1 which 

verifies 

M 

U < c£ n2 ' i “ ^ 2 - < 2 - 2 ) 

»=1 

Here C is a positive constant independent of N and 0. Defining = rot 0 at we 
have that w^r € and 

||w — Wjvll! < ||0 — 0jy|| 2 ^ l[0*lL;.i — CT. N 2 **' 

i=l »=1 

for m > 2. This gives (2.1) for Oi > 1. I 

Theorem 2.2. Let u and ujy be the solutions of problems (1.2) and (1.4) respec- 
tively. If u l € (H*' (Cli )) 2 , > 1 and f € (H 71 * (fi»)) 2 , /a* > 1, for i = 1 then 

there exists a positive constant C independent of N, u, and f such that 

||u — ujy ||i < C { N l 9% |iu t || <y . i + iV 1 (2*3) 

i=l 
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Proof. From a general approximation result (see e.g. [3, Corollary IE.3]) we have 


1 


<c{ 


sup 

*n€Vs 


=0 {ll u-w "H 1+ 

((Vw//, Vzjv)) - ((Vwjv, Vz N ))u 


inf 

w w €V , Wl V*w w 


(2.4) 


l|Vz*|| 

((f,Zj/)) - ((f,Zjy))^ 


+ 


■i»€?w IIVzatII 


}■ 


where C is a positive constant related to the constant Ci appearing in (0.5) . Choosing 
w u = IlJ^u the first term at the right hand side of (2.4) can be estimated using 
(2.1). Moreover, since the Legendre Gauss Lobatto quadrature formula (see (0.1)) is 
exact on i 3 2 JV-i(Gi), the second term vanishes. The last term in (2.4) is estimated 
by using known results on the interpolation error (see [8, section 9.4.3]). I 

Our aim is now to prove an error estimate for the pressure. We remark that if 
Mfj is given by (1.12) then pu vanishes at the comer points ao, bo, a m, b\{. On the 
other hand the solution p of the continuous problem do not vanishes at these points. 
Hence if Mu is defined by (1.12) the spectral accuracy cannot be achieved. Using 
an idea introduced in [6] we now define a new pressure space Mu : The following is 
a direct consequence of theorem 1.2. 

Proposition 2.1. Let Mu be a subspace of Il^x ^V(G») such that 


Qn,x := {g € { IJ P[a^j(U»)}/R : 9* = 9* +l at ai and &», i = 1,...,AT - 1 j C Mu 

»= l 

(2.5) 

where A € (0, 1) and [A N] is the integer part of A N. Moreover suppose that there 
exist an isomorphism Gu : Mu — ► Mu such that for any q € Mu 


((9> V • v))at = {[G^q, V -v))// Vv € V N (2.6) 

and 

\\G N q\\ < C\\q\\ (2.7) 

with a constant C independent of q and N. Then the inf-sup condition (1.13) holds 

with Mu substituted by Mu. I 

A new discrete problem is obtained by using Mu instead of Mu in (1.4). From 
the proposition 2.1 it follows that such problem has a unique solution (u^ ,pjv) and 
that ujv = u/v. In other words the discrete velocity do not depends on the choice 
of the pressure space provided the assumptions (2.5)-(2.7) hold. We now give an 
example of pressure space. 

Given an integer i, 1 < t < M, and a point P € E* N , the characteristic function 
Xp (defined in (0.6)) can be written in the form 

N N 

Xp = ^2 ^ a k,iLk{x)Li{y)- (2.8) 

k=0 1=0 
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We define Xp by 


a- t Y2 a ULi( x ) L i{y)> ( 2 * 9 ) 

fc=[Atf] + l /=[AAT ]+1 

and we consider the space: 


M 

1=1 

(jSxi.Jw.l = («‘>XfcJw, l = (« W ,x " ,)nm = («".xJ[ ( )w.M =0, 
i‘(h) = ? i+1 (*«i), >' = 1 M - 1, 

M M 

1)^,1 = 0, y^(g*,^(y))j/,i = Q, 

t-i »=i 

{qi'L^ix))^ = 0, {q\U N {x)L N {y)) N ,i = 0, * = l,...,Af}. (2.10) 

Recalling (1.12) we remark that the difference between Ms and Ms consists_in 
the first four conditions appearing in (2.10). From (2.8)-(2.10) it is clew that Ms 
verifies (2.5). We now prove the following. 

Proposition 2.2 If Ms is given by (2.10) then there exists an isomorphism Gs 
such that (2.6) and (2.7) hold. 

Proof. Let us denote by Xfc> k — 1, ...,4 the four elements of fli^i Pn[^*) given 
in (1.6) (1.7) and by Xfc, k = 1 ,..., 4 the corresponding elements with Xp substituted 
by Xp (see(2.9)). We define Xk, k = 1,...,4 as the projection, for the inner product 
(M)*, of Xfc on the subspace of Pjv(Hj) which is orthogonal to the elements 
given in (1.8)-(1.10). 

For any q € Ms we define 


4 

Gsq = <1 ~Y1 

fc=i 


((g,Xfc))jr 


(2.11) 


From the definition of Xfci k = 1, ..., 4 we have that Gsq is orthogonal to the 
elements (1.8)-(1.10) for any q € Ms- Moreover using the expansion of Xfc, Xfc, Xfc, 
1 < fc < 4, in Legendre series (see e.g. [6, Prop. V.3]) one can show that 


((Xfc,xO)w = and ||Xfc||WIIXfc|_k < C k,l = 1 ,..., 4 , (2.12) 

where 6k, i is the Kronecher symbol. Using (2.11) it follows that Gsq is orthogonal 
to X/, / = 1,...,4, hence Gsq € Ms for any q € Ms- 

Since Xfc € Zs,m k = 1,...,4 we deduce from proposition 1.2 that (2.6) holds. 
To prove (2.7) we remark that using the Cauchy-Schwarz inequality and (2.12) we 
have 

BGirtlU < MU (i + E j||j|£) < CM*. (2.12) 
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Hence (2.7) follows from (0.5). 

We can now give a convergence estimate for the pressure. 

Theorem 2.S. Let (un,Pn) be the solution of problem (1.4) with Mn sub- 
stituted by Mn- If (u,p) is the solution of problem (1.2) and u‘ € (H 9i (H»)) 2 , 
p e H‘ r< “ 1 (n») n S l ( n), > 2, f € (H^n*)) 2 , tii > 1, 1 < * < M, then we have 

f)Ar j -'<{iiu i iu,i+iiP i iu- 1 , i }+f;jv 2 -»<iiriu, i }, (2.i3) 

»= 1 


Up - pwIIo < C'j 


where C is a positive constant independent of N, p, u and f. 

Proof. From a classical approximation result (see e.g. [3, corollary 2.3]) and 
using theorems 1.2, 2.1 and proposition 2.1 we obtain 


IIp-pk||o < civ{ 


Af 

inf l|p-? W || 1 , + £{iV l -'‘||u i |U ( . j +W‘-‘‘<||f|| 

<7n6A Is 


«,<}}• ( 2 - 14 ) 


Recalling (2.5) and an approximation result due to Bemardi and Maday (see [5 thm. 
IH.3]), we obtain 


M 

inf ||p-«w||o < |nf II? - «rl|o < C^W^'llp'k.i. ( 2 - 15 ) 

<jn€M s 1 n€Qn,x f=1 

hence (2.13) follows from (2.14) and (2.15). I 

Remark 2.1. Let (vln,Pn) be the solution of problem (1.4), and let (u n,Pn) be 
the solution of the same problem with Mn replaced by Mn- Recalling (2.11) (1.4) 
and (2.6) it is readly seen that 

PN = GnPn- (2.16) 

Hence to obtain a spectral approximation of the pressure one can either sojve (1.4) 
and compute pn from (2.16), or solve directly (1.4) with Mn replaced by Mn- I 
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3. Plurirectangular Domains. 

In this section we consider a plurirectangle 0 partitioned in rectangles 0<, i = 
1, ...,M, where M is a fixed integer. We assume that the decomposition of 0 is such 
that the intersection of two adiacent subdomains is either a point or an entire edge. 

The finite dimensional space V n is defined again by (1.3) and the discrete prob- 
lem is given by (1.4). The crucial difference with the case of aligned subdomains lies 
in the choice of the pressure space M/y. Let us first consider, as an example, the 
domain of figure 3.1. 



figure 3.1 


The points a», t = 1,...,6, are the non-reentramt corners of 50; , * = 1,...,6, 

are the points of di 1 which are the intersection of two subdomains, ci and C2 are the 
reentrant comers of <90 and d t is a so called “cross point”, i.e. the intersection of 
four subdomains. 

Following the lines of section one we define Mu as the subspace of nj=i -Pw(Ot) 
whose elements are orthogonal to (see (1.5) and (i.ll)). In this case there are 
30 independent elements of 7, and the space Mu is given by 

r 

M n = (g € JJ P^(Oi) : ? = 0ma<, t = l,...,6, (3.1) 

»=i 

q' = q *+ l in b { , i = 1 ,..., 6 , 

q l +q 5 = q 4 in a, q 5 + q 7 = q Q in c 2 , 

q 1 + q 3 = q 2 + q 4 in d 1, 

7 

y~!(g*» i)jv,» = 0, 

i— 1 

((?%■£/* ( x )£j\r(y))-tf,* = 0 * = i,...,7, 

2 6 

{< 1 7 > l n{ x ))n , 7 = 0 , y>,LWz)W = 0 , VXt'LM*))," = 0 , 

»=i i= 3 
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7 

y^(g*>-Ew(v))^.» = o> (9 5 >^(y))/^,5 = 0) (9 1 >^(y))^,i + (9 4 >^(y))iv,4 = o» 
*=6 

5^(g*> £jv(y))w,» =o}. 

i=2 

The first 15 conditions are orthogonality relations, as a matter of fact they can be 
written as follows. 

(?*> Xa,-)-Ni» = J = 1,2,3, (q 6 ~ »Xa#)^.7 (? >Xa*)-f/,7 

(y'jXij)/^,! = (g* +l > Xbj 1 " 1 )^,!? * = 

(ySxcJ^a + (? 5 ,XcJtf,5 = (g 4 ,x2Jw,4, ( 3 -2) 

(? 5 ,XcJ*,5 + (q 7 ,xl 3 )N, 7 = (q 6 ,Xc 7 )x,6, 

(gSxaJwa + (g 3 .XdJjv, s = (g 2 ,XdJ*,2 + (g 4 ,XdJjM- 

These conditions will be called “physical conditions” as opposite to the last 15 con- 
ditions appearing in (3.1) which will be called “spectral”. 

We now suppose that 0 is any plurirectangle partitioned in M rectangles 0,, 

* = 1, ..., M. Let Mn be defined by (1.11), in the following we describe the conditions 
which characterize the elements of Mu . 

A physical condition is associated to every node of the decomposition. Precisely 
if P is a non-reentrant corner of 30 (as a,- , t = 1, ..., 6 in fig. 3.1) then we impose that 
any q 6 Mu vanishes at P. If P is a reentrant comer (see fig. 3.2) then we impose 
5*(P) + ,“{P) = ^(P). 


p 



0 n 



°k 








figure 3.2 

If a node P belongs to 30 and it is not a corner (as 6,-, i = 1, ...,6 in fig. 3.1) we 
impose the continuity of q at P. Finally if P is a cross point (see fig.3.3) we impose 
the “continuity” condition g m (P) + ^(P) = g*(P) + g n (P). 

Let us now describe the spectral conditions appearing in the definition of Mu • As 
in the continuous case we impose that any q 6 Mu is orthogonal to the constants. The 
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P 

c 

G 




0, 



figure 3.3 


orthogonality of q to L i N {x)L i N {y), i = gives Af additional conditions. Let 

us now call “horizontal strip” any maximal set of horizontally aligned subdomains. 
For instance in figure 3.1 these strips are {^ 7 }, {fii,ft 2 }and{fL, * = 3,..., 6 }. If a 
horizontal strip is given by {fL, / € A} where A is an index set, then we impose 
that q is orthogonal to the function which equals L'n (y) on flj, / € A, and vanishes 
on the other subdomains. Similarly if {f2 fc , k E K} is a “vertical strip” (these 
strips are {n 6 ,n 7 }, {n 5 }, {n^rUjand^.Hs} in fig.3.1) then we impose that q is 
orthogonal to the function which equals L 1 ^ (x) on fl*, k E K, and vanishes on the 
other subdomains. 

Using this definition of the pressure space Mn, the “inf-sup condition” (1.13) 
holds also in the case of a plurirectangle (see remark A.l). In particular this im- 
plies that the discrete problem is well posed. Moreover, it is readly seen that if 
w N G VjythenV • Vn E Mjv, so that the discrete solution is divergence-free for any 
value of the discretization parameter. 

We now give the collocation form of problem (1.4) for a decomposition of a 
plurirectangle fl. 

Let (ujv»Pjv) be the solution of (1.4), we have 

-Auj,(£) + Vp-„(£) = f(£), V Hi, n n ( , = 1 M. (3.3) 

Denoting by the common side of two subdomains fl*, flj , (not including the 
end-points) we have 

< 3 - 4 > 

= - { - Au J, + VpJ, -f‘}(£)^ - { - Au'„ + Vp'„ -f'}(£)^, V{ € I*,,. 

Let £ be a cross-point and consider, for simplicity the notation of figure 3.4. 
Setting 

DH = E -pW’*}(«) 

k=l 
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figure 3.4 


the cross-point condition is 


£ pWHO = -£{-au' n + v„v,-r}({)<4, 

*=1 »— 1 


V cross — point £. (3.5) 


Finally Um verifies the incompressibility condition, i.e. 


V • u' N =0 in Z. % N , t = 1, (3.6) 

Actually, it is sufficient to impose (3.6) in a subset of {B^, * = 1 , ...,M}. Precisely, 
since um = 0 on 30 then V • um automatically vanishes at any non-reentrant corner 
of 30. If P is a node on 30 which is not a comer (e.g. 6i in figure 4.1) it is 
sufficient to impose only one condition at P (e.g. for fig. 4.1 V • u]y(&i) = 0 implies 
V • u^(&i) = 0 since V • um € Mm). Similarly if P is given in figure 3.2 or 3.3 it 
is sufficient to impose V • = V • = 0 and V • = V • — V • = 0 

at P respectively. It is also not necessary to impose (3.6) at some collocation points 
internal to each subdomain (this follows from the fact that V • Um € Mm and hence 
it verifies the spectral conditions). Precisely one can eliminate in (3.6) one point for 
each subdomain, one point for each horizontal strip, one point for each vertical strip 
and one point choosen in an arbitrary subdomain. The points which are eliminated 
must verify some hypothesis of the kind of (1.8), (1.9); if in every subdomain these 
points are not aligned then these hypothesis hold. Summarizing, (3.6) reduces to 

V-u 1 m = 0 in Bm/Am.j w, * = 1,...,M 

where Am , m is a set of n collocation points and n is the codimension of Mm in 

The velocity error estimate (2.3) holds also in the case of a plurirectangle. The 
proof works exactly in the same way. 

As in the case of the aligned subdomains to obtain an error estimate for the 
pressure one has to modify the space Mm introducing a new space Mm such that 
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the thesis of proposition 2.1 holds. Precisely the physical conditions related to the 
corners of dCl have to be removed. In the case of a non-reentrant corner P€ 517* 
the condition 9 *(P) = 0 has to be substituted by {q*,Xp)N,i = 0 (as in the case 
of the aligned subdomains, see (2.8)-(2.10)). On the other hand if P is a reentrant 
corner (see fig. 3.2) the condition q k + q n = g* at P has to be substituted by 
(g*>Xp)i^fc + (9 n >Xp)jv,ri = (q 1 , Xp)n,i- If the other conditions of M N are not changed 
we obtain a new space Mn which verifies (2.6), (2. 7) and: 

M 

€ { U P[xjv](iit}/I^ : Q is continuous at the verteces of the subdomains j C Mfj- 

Then the result of proposition 2.1 holds and the estimate (2.13) can be obtained with 
exactly the same proof technique. 
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Appendix 

In this section we prove proposition 1.1 and theorem 1.2. We first consider 
two aligned subdomains, the proof in the general case will follow from an induction 
argument. Let us precise some notation. The outward normal unit vector to dQ 
is denoted by n. For i = 1,2 we denote by r* the counterclockwise unit tangent 
vector to dUi and by n* the outward normal unit vector to The notation 

rj, n*, i,j = 1,2, are introduced in figure A.l. 



fig. A.l The decomposition of fl with M = 2. 


Recalling the definition (0.6) we set for any n > 1 : 


Jn 


L>n+ 1 ~ Ln-l 

2n + 1 

—■ ZJizl . Ss±l ~ L n~ . i fort = 1,2 

2 2n + 1 


where {x^} are the abscisses of the points {oi} for t = 0,1,2. The polynomials J n 
and J\ are the primitive of L n and L\ respectively, i.e. 


j; = L n , J„(-l) = Jn(l)=0 

(4)' = Ln> = «£(*) = 0, fort = 1,2, 


Moreover we have for any n > 1 and t = 1,2 : 


I T* ||2 __£l_ 

I K Ho,.* - 2 n + 1 


( 2 » + 1 ) 


3 < IK Wb 


(2n + 1)' 


(A.2) 
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where C \ , C 2 and C$ are positive constants independent of n. 

Let us denote by Zff, 2 the subspace of x Pw(U 2 ) spanned by the 12 

elements _ ... 

(xi„,0), (xj„, 0), (0 ,xi), (0,x2,), (A-3) 

(x^.-xi.), (xi,.-xl,). (A-4) 

(1,1), (Lr,(y),L K (y)) (A.5) 

(L# (x),0), (0,L%(x)), (L l s [x),L N (y), 0), (0,L|,(i)£„(y)). (A.6) 

These are precisely the spurious modes given in (l.6)-(1.10) for two subdomains. 
Using (0.2) it is readly seen that Z ^,2 is contained in the space Zs , 2 defined in (1.5). 
In the sequel we prove that, actually, Zn , 2 and Zn ,2 coincide. 

Let q € Pv(fti) x P N {U 2 ) be orthogonal, for the inner product (0.4), to Z N , 2 
(i.e. to the elements (A.3)-(A.6)). For i = 1,2, q i can be written in the form 


N - 1 


N - 1 


(A.7) 


(A.8) 

(A.9) 


q* = ^ a ki^k ® Li + ^ Pk^'k ® «^r-i + ® P/- 

0<fc,/<iV-l fc=o <=l 

Let us define two continuous functions A and B: [xo, x 2 ] — ► R such that 

A and B are linear in [xi_i, x<] , « = 1, 2 

A(x 0 ) = P(xo) 

and 

Vz €]xi_ 1 , Xi [ Al(x) = <4 s, ( x ) = 4 for * = h 2. 

We define the function v = (t>i,u 2 ) € {Pn{^i) x Pn^)) 2 as follows: 

N- 1 

V* = 4/4 44 ® J N-1 + A* ® Lo + B i ® Jat-1, 

0</<fc<l'/-l,fc+/*0 fc=l 

(A.11) 

JST— I 

«2 = ® *7/ + yi '//4r~i ® *^/ ^° r * = 

0<k<l<N-l /=1 

Lemma A.l . Let g € P//(Hi) x P^^) be orthogonal to ijv,2 for the inner 
product ((.,•))* defined in (0.4). If q is written in the form (A.7) and v is given by 
(A.ll) then we have 


v* = q x for i = 1, 2, 

(A. 12 ) 

v • n = 0 on dCl 

(A.13) 

4 = -v 2 • n 2 on Ti. 

(A. 14) 
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Proof. From (A.l) (A.7) and (A.ll) we obtain (A.12). Using (A.l) we deduce 
that v • n = 0 on the two segments [a<j, a?] ^d [&o> 62]" Moreover we remark that 
the orthogonality of q to the elements given in (A.5) reads as follows 

2 2 

t=i i=i 

From (A.8)-(A.10) we deduce 


A(x 2 ) = B(x2) = 0. 


(A.15) 


Then the condition v • n = 0 on the two segments [<zo» ^0] an d [a2, 62] follows from 
(A.1)(A.9) and (A.15). Finally since A and B are continuous we deduce (A. 14) from 
(A.l) and (A.ll). 

Lemma A.2. In the same hypothesis of lemma A.l the following estimate holds 


Dk lit, < 


(A. 16) 


i=i 


i=l 


for a positive constant C independent of q, v and N. 

Proof. Using (A.2) and (A.7) a direct computation gives, for i = 1,2: 

ll9‘llo,i ££((*< E (“il) a (2t + 1)(2J + 1) + 


(A.17) 


0<kJ<N-l 


N 3 


, / \3 N-l , 

1 , (*i ~ ^i-l) V' ( ~i\ 2 1 

N 3 


2k + 1 


+ " E (Tf^a + l)’ 


fc=0 " 1 = 1 

where C is a positive constant independent of q , N and i. It is readily seen that 


^ c'lklo.o * - 1 > 2 * 


(A.i8) 


r\ I 

Let us now estimate j| ||o,» * From (A.l) and (A.ll) we obtain 


dv\ 

dy 


N - 1 


^2 d kl J' k ®L\+Y, Pi4 ® L N- 1 + B i ® Ljv-i- [A. 19) 


0<l<k<N-l, k+l * 0 fc=l 

We now estimate the three terms appearing in the right hand side of (A. 19). 
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1) Recalling the inequality {L\) 2 dx < C(2l + 1), we get for * = 1, 2 


E “44 ® < c [ e ( 2i + 1) ( r ( £ °u4) 2 d*) 

0<l<k<N - 1 1=0 Jx *-i k=l, k+l* 0 


1/2 


< 


(=0 

N-i 


k=l t 


s c[ £ (a * 

Using the Cauchy-Schwarz inequality and (A. 17) we deduce 


ii ^2 a ki^k ® -^<iio,t ^ 

0</<*<^-l 

N - 1 tf-1 

<C(E2' + i)(E(“«) 


(A.20) 


1=0 


k= 0 


,i \2 (j* J z-l) 3 

(2*+ 1)(2/ - 1) 


) < ClV 2 || g i || 2 i 


2) From (A.2) and (A.17) we obtain for * = 1,2 


(A.21) 


II E A4 ® < f E (^) 2 ^ -- + i) 1 i ;i g^ll^llv- 

3) Recalling (A.9) and using the Poincare inequality, (A.10) and (A.17) we deduce 
2 2 

]T || B* ® L^-ilH.i < f / ’ B 2 (x)dx < ^ - x,_0(/?‘) 2 .< (A. 22) 

*=1 Jx <> V i=l 

< CiV 3 2 II, <11^. 

From (A.19)-(A.22) we obtain 


Ell$fll2, i <cjv 2 Elk'll^. 

1=1 y i= 1 


Using the same technique it is possible to estimate Y^i=i Iloilo,* > hence, recalling 
(A.18) we obtain (A. 16). | 

Lemma A. 3. In the same hypothesis of lemma A.l the following estimate holds 

2 2 

J2 H v< ' r< IU»(ao.) 1 ^ c J2 ( A - 23 ) 


1=1 


i= 1 
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for a positive constant independent of q, v and N. 

Proof. From (A.l) and (A. 11) we obtain for i = 1,2 

V*€ v* • rl(x, l) = — E ««•**(*) + «**(*)• (A.24) 

0 <l<k<N — 1 k +&0 

Using (A. 2), the Cauchy- Schwarz inequality and (A.17) we deduce 

N-i _ k 


r ( E <,4(*)) 2 dx < c £ [(£ 4,) 2 %r 

Jxi ~ l 0<l<k<N-l, k+l* 0 k=0 1=0 \ K 


,i 

+ 1) 3 


N-i 


< C{xi-Xi-if E [(E 


a 


kl I 


k 

N-l k 


) 2 ~r~ — 1 < 
' 2k + 1 J “ 


H ( 2 ^ + 1 ) 1/2 ( 2 ^ + 1 ) l/2; 2Jfc + 

(4/) 2 k + 1 


< C(x, - *<-i) 3 E E (^ i)( 2 / + i )2fc + +1 * wwb- 

The second term at the right hand side of (A.24) can be estimated using the 
technique of the proof of (A.22), so that we obtain 

2 2 

E« vi - r ilU.(™o<cElk‘llS..- 


1=1 1=1 

Working in the same lines to estimate ^Li ||v i • j = 1,2,3. we obtain 

(A.24). , 

Theorem A.l. Let q € Pjy(Hi) x Pn (f^) be orthogonal to Zs, 2 . There exists 
w € Vfj such that 

V ■ w‘ = q' ,for * = 1,2 


and 

E iiw‘iii,i <cN 2 f2 iisio.i , 

1=1 1=1 

for a positive constant C independent of w, q and N. 
Proof. Let and 0# be two functions such that 

*j. &N € Pjv([— 1 , 1 ]) 

«J(±1) = 0, «;(±1) = 0 

W)'(-1) = 0, (^)'(-i) = 1 

(9J)'(1) = 1, «S<(1) = 0 


(A. 25) 
(A.26) 


(A.27) 
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and 


(A.28) 


lWM-i. 1 ) + ll®JrlU*(-M) - jv3- 

One can for instance choose 0^(y) = — ^ ^jv(y) = ~ J/)- 

(it can be proved that it is not possible to find 0% and Ojj verifying (A.27) and (A.28) 
with on the right hand side a power of N greater than 3). 

Let q and v be given by (A.7) and (A.ll) respectively. We define the functions 

/ and g as follows: 

/(y) = v 1 * r i (&i)0jv(y) -yl * r a( a i)*ir(tf) 

(A. 29) 

»(») = -^(MWv) - f“ v s [-l.i]. 

For * = 1, 2 we now define the functions <j> x and t p* on dfl» as follows: 

rpi = 0 on dCli Pi dfl 


Vy e [—1,1] ^(xi ,y)=/(y) 

_ yi . T » Qn n dn 

Vy € [-1, 1], = g(y) and i> 2 {x t ,y) = v l • rj + v 2 • if - g{y). 

We now show that there exists $ 6 Pn{^i) x -FW(^ 2) suc h that 

= <f? on ddi 


^7 = on an t - 


and 


5^11**111 ,< - ^^(^ll^llL^dOi) +^ v ’ 2 H^*IU a («n<))* (A. 32) 

»=i »=i 

The existence of this function $ is a consequence of a trace theorem due to Bemardi 
and Maday provided some compatility conditions hold (see [4, thm 3.1]). To simplify 
the notation let us fix the index t and set Ai = Oi_i, A2 = A 3 = i>i, A4 = &»- 1 
(see fig. A.l). We denote by <j>) (resp. rjj)) the restriction of <f>' (resp. tp') to [Ay_i, Ay]. 
The compatibility condition are the following 


4(A j )=4 + 1 U j ), 1<J<4 

07 i 

¥&(*/)--$&). 1<3<* 


di>) 


dr) 






1 < j < 4. 


(A. 33) 
(A.34) 

(A.35) 

(A.36) 
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Thanks to (A.27) (A.29) and (A.30) the conditions (A.33) are satisfied. 
For j = 1 the condition (A. 34) reads as follows: 

v 1 •T’ 2 1 ( a o) = 0 

v 2 *r|(a 1 ) = -/ / (-l). 

Due to (A.27) and (A.29) these relations can be written in the form 

v 1 • nj(oo) =0 


v 3 -n^fli) = -Vi -n^aj) 

hence they follow from (A.13) and (A. 14). Using (A.13) (A.27) and (A.29) and 
working in the same lines the conditions (A.34) for j = 1,...,4 and (A.35) can be 
easily verified. For j = 1 the compatibility condition (A.36) reads as follows: 


dv l • Tj 
dT l 


(oo) = 


dv 1 • rj 


M 


dv 1 • ri 
dr 2 


(®i) + 


dv 2 ■ T 2 

dr 2 


K) + s'(-i) 


dv 2 ■ r 3 


(«i). 


The first condition can be written in the form V • v 1 (a 0 ) = 0. It follows from (A.12) 
and the orthogonality of q 1 to xi 0 • Recalling (A.29) and (A.30) the second condition 
is V • v x (ai) = V • v 2 (ai). This is a consequence of (A.12) and the orthogonality 
of q to (xi^-XaJ- Working in the same lines the relations (A.36) for j = 2,3,4 
can be proved. They follow from (A.12) and the orthogonality of q to the elements 
(°>X2,)> (X6 o »0), (0, X 2 J and (x^.-X^)- 

It follows from (A.33)-(A.36) that there exists $ € Pjv(Dx) x Pn(C12) such that 
(A.31) and (A.32) hold. We set 


w = v + rot$ , (A.37) 

where v is given by (A.ll). From lemma A.l and (A.30) (A.31) we have that w € V N 
and V • w = q. To prove (A.26) note that 

H H w< Hm ^ c'EOK’llli + \\Hb)- (a.38) 

*=1 »=1 

Recalling (A.32) we have to estimate the terms ||^ i |lL J (an < ) J and 

Yli = i llV'NI.L^dn,)’* Thanks to (A. 28) (A.29) and (A.30) we have 

^||^||L»(ao,) < Tr6(l(v 1 T‘)(6 1 )| 2 + |(v 1 -r 2 1 )(a 1 )j 2 ) (A.39) 

»=i 
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and 


£ < ^d^r i (**)l a +l 2 ^T 2 (««)l*)+£ h'-W* (>•«) 

From (A.39) (A.l) (A.10) (A.ll) and (A.17) we get 


, dv 1 • r. 1 


< Jjs\ Ai ( x i )\ 2 ^ jye l^ool 2 ^ ^6lk l |lo,i 

t=i 

From (A.l) (A.ll) and (A.40) we obtain 

X^II^IU s (3n<) < S&( S l-Llf + EK^M-Oi)* 

i=l 0</<Jk<AT-l i=l 


(A.41) 


Using the Cauchy-Schwarz inequality, (A.17) and lemma A.3 we deduce 

£||^IIl.,™,) < ^ os| E k _. W‘l 2 (2i + i) 1 (a + i) + {AA2) 

+ c"£l|« i ll?, i <c£||, i ||?, j . 

1=1 »=i 

Finally from (A.37), (A. 16), (A.38), (A.32), (A.41) and (A.42) we obtain (A.26). I 

We are now in a position to generalize the result of theorem A.l to the case of M 
aligned subdomains. To this end let us denote by Zn,m the subspace of n»=i -P/v(U») 
spanned by the 4 M + 4 spurious modes given in (1.6)-(1.10). The other notation axe 
those of section one. 

Theorem A. 2. Let q e II^i -Pjv (^i) be orthogonal, for the inner product ((•, -)) N 
defined in (0.4), to any element of Z N , M . Then there exists wgV w such that 


V • w 1 = q* for i = 1, . . . , M 


(A.43) 


and 


M 


M 


£l|w i ||?, i <cjv 2 ^||, i ||§, i , 


(A.44) 


i= l 


i= 1 


where C is a positive constant independent of w, q and N. 

Proof. We use an induction argument on M. Theorem A.l gives the case M = 2. 
We suppose now that the theorem holds for M — 1 (Af < 3) aligned subdomains and 


24 


we consider the case of figure 1.1. The function q 6 FI^i -Pv (fi*) can be decomposed 
as follows: 

g = g + ? (A.45) 

where 

g 1 = q 1 , g* = o for i = 3, M (.4.46), 

and g 2 has the following properties 

f(ai) = q l (ai), q 2 {h) = g 1 ^), g 2 ^) = g 2 (& 2 ) = 0, (A.47) 

{q X iI&Lo)s %x + {q 2 ,L 2 Lo) N 2 = 0, (g 1 ,ijv(y))j\r,i + (g 2 »£/v(y))tf ,2 = 0, 

(4.48) 

{q 2 ,L 2 N {x)) N2 = (q 2 y L%{x)L N {y)) N2 = 0, 
there exists a positive constant C independent of N such that 

llfllw < C\W l || M . (A.49) 

The polynomial f 2 is then defined by ^ = g 2 — g 2 . Let us assume for a moment that 
such a decomposition exists. From (A.45)-(A.48) it follows that (g 1 , g 2 ) is orthogonal 
to Ztf '2 and (g , g 3 , . . . , q M ) verifies the hypothesis of the theorem in the M — 1 
subdomains > • ■ • , fW- Hence from theorem A.l and the induction hypothesis there 
exists W £ nfssi Pw(^») 2 and ^ € \\iL 2 Pw(^») 2 with w = 0 on 3(fli U = 0 

on d(u££ 2 fi»), w and W continuous at the interfaces, such that 


V-w< =g\ * = 1,2, 

II 

> 

g*, t = 2, . . . , M 

(A.50) 

and 




E K'liw < cjv 2 E nfiio.i 
1=1 1=1 

M 

. Ell*l 

t=2 

IwSCJV^IlfllL. 

i=2 

(A. 51) 

We now set 

f w 1 m 
w = < w 2 + W 2 
in 

\ 

Hi 

in fl 2 
Hi for t 

= 3, . . . , jlf. 

(A. 52) 


Clearly w E Vjf , furthmore from (A.45) and (A.50) we have that V w‘ = g* for * = 
1 , . . . , M. 

From (A. 52) (A. 51) and (A.46) we obtain 


M MM 

E iK'iil- * c { E lie'll ?.< + E ii*'iiw) £ 

»=i »=i %= 2 

M 

<CN\ E ll«*llw + + Ilf Ilia) 

>— 1 2 
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hence using (A.45) and (A.49) we get 


E iKllii < CN* 
1=1 


( E Il? < llo,i + ll9 1 |lo.i + 

i=l, ijt 2 



which gives ( A.44) . To conclude the proof we exhibit a polynomial q 2 which verifies 
(A.47) (A.48) and (A.49). We suppose, for simplicity, N even and we set: 

f =K) 2 X?, 9 l (!’ 1 ) + K) 2 X«y(“i) + <*.Xg ® io + a 2 Ll ® X„+ (X. 53 ) 
+ 03 Ljy ® Xo + o . Lifj ® L w + 03X3 ® X3 


where a*, A: = 1, ... ,5 , are constants. 

From (0.7) it follows that (A.47) is verified provided as = — ^ fc=1 otk • 

An easy computation shows that ak, k = 1, — ,4, can be found so that (A.48) 
holds, moreover we have 

|«fc| < (IkMkl + (Wjv) 2 |tf l (&l)|||xi> x ||o,2 + ( w Ar) 2 | 9 1 ( a l)IIIXail|o, 2 )- 
Hence from (A.53) we deduce 

||g 2 ||o,2 < C'(lk 1 ||o,i + ( w Ar) 2 (l9 1 (^i)lllx 2 l ||o,2 + |9 1 ( a i)|||x 2 l ||o,2))- 


Using (0.2) and (0.3) it is readly seen that the two last terms can be bounded by 
H^llo, 1 hence (A.49) holds. I 


Proof of Proposition 1.1 and Theorem 1.2. 

From theorem A.2 and (0.5) we have that there exists a positive constant C 
. independent of N such that 


Vq e Znm 3 w € V N 


((<?) Y . ' w ))jy r Ikllo > C_ 
IklUIIVwlU - ||Vw||o - N' 


hence 


inf 


(k,V-v))w ' g 
ll«IUI|Vv|| w - N 


{A. 54) 


From (A.54) we deduce that Z £ M do not contain any spurious mode for the pressure. 
Recalling (1.5) we then have Zjj M D Zn, m = {0}. This implies that 

dim Zif.M < dimZjv.A/* Since Z^,m D Znm we deduce that Zn,m — Zn,Mi hence 
proposition 1.1 is proved and theorem 1.2 follows from (A.54) and (1.11). I 
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Remark A.l (More complex decompositions). 

The proof technique of theorem A. 2 can be generalized to cover more complex 
decompositions. Consider for instance the L-shaped domain of figure A.2.a. 



figure A.2.a figure A.2.b 

If q £ Mu (the pressure space described in section 3) then it can be written 
in the form (A.45) with q 1 = q l , q 2 = q 3 and q 2 such that (q 1 ,q 2 ) belongs to 
Zjj 2 - ft turns out that theorem A.l can be applied in fli U Q2 and in U fl3 
and that theorem A. 2 holds for this decomposition. Analogous consideration holds 
for the figure A.2.b. In this case q € M ^ can be written in the form q = q + q 
with q l = 0, q^ — q^i <q^ = 9 1 1 qt — 0 and q 1 , q 1 for t = 2,3 such that one can 
apply theorem A.2 on the L-shaped domains Hi U fl2 U fi3 and fi2 U 03 U Q4 . For a 
general partition of the domain of the kind of that of figure 4.1 one can always find a 
decomposition of q as in the previous cases. Hence proposition 1.1 and theorem 1.2 
can be proved also in this case. I 
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